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ABSTRACT 


This  paper  presents  the  development  and  evaluation  of  two 
algorithms  for  the  determination  of  wind  compensating 
launcher  settings  for  unguided  rockets.  The  algorithms 
are  designed  for  use  in  an  iterative  trajectory  simulation 
process  and  can  be  used  to  determine  launch  angles  which 
will  yield  one  of  many  possible  trajectory  objectives. 

One  algorithm  provides  the  first  estimate  as  a  function 
of  the  ballistic  wind;  the  other  controls  the  iteration. 

Four  test  cases  are  presented  to  evaluate  the  technique. 
Sixteen  different  wind  profiles  were  applied  to  four 
separate  unguided  rockets.  Two  different  trajectory 
objectives  were  considered  -  nominal  burnout  attitude 
and  nominal  impact.  The  tests  show  that  the  process 
converges  rapidly,  usually  requiring  only  one  or  two 
iters tions. 
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INTRODUCTION 


The  trajectory  of  an  unguided  rocket  or  projectile  will  be  affected 
by  the  wind  encountered  during  flight.  This  pertubation  of  the  tra¬ 
jectory  is  commonly  calxed  wind  effect.  If  one  desires  to  obtain  a 
predetermined  value  of  a  given  trajectory  objective  (such  as  impact 
of  a  given  booster,  orientation  of  the  velocity  vector  at  burnout, 
etc.),  it  is  necessary  to  determine  a  launcher  setting  which  will 
compensate  for  the  wind  effect.  The  sophistication  of  the  algorithms 
to  determine  these  launcher  settings  depends  not  only  on  the  specified 
trajectory  objective  but  also  on  the  rocket  configuration  and  the 
nominal  trajectory. 

A  first-order  approximation  of  the  wind  effect  may  be  obtained  by  a 
technique  called  ’’wind  weighting."  This  technique  is  discussed  in 
Appendix  A.  Several  procedures  [Lewis  (1949),  James  and  Harris  (1960), 
Hennigh  (1964),  Duncan  and  Engebos  (1966a,  1966b)]  based  on  the  wind 
weighting  technique  have  been  developed  which  give  first-order  ap¬ 
proximations  for  the  launcher  settings  required  to  compensate  for  the 
wind  effect. 

In  the  previous  developments,  the  model  was  usually  restricted  to 
considering  either  a  specific  trajectory  objective  or  a  specific 
type  of  trajectory  or  both.  This  paper  presents  a  procedure  which 
allows  for  a  much  wider  range  of  applicability. 

The  model  presented  herein  has  a  twofold  capability.  An  approxima¬ 
tion  for  launcher  settings  based  upon  the  wind  weighting  technique 
is  developed.  For  missile  projects  which  require  more  accurate  com¬ 
putations,  an  iterative  procedure  is  developed.  The  algorithms  for 
the  iteration  converge  rapidly;  the  wind-weighting-based  model  is 
used  as  a  first  approximation.  The  iterative  procedure  can  also 
be  used  to  refine  the  numerical  values  for  the  wind-weighting-based 
model.  The  use  of  the  complete  model  as  an  operational  tool  re¬ 
quires  a  real-time  computational  capability  (Duncan  and  Rachele, 

1967). 


COORDINATE  SYSTEMS  AND  TRANSFORMATIONS 


In  this  development  it  is  assumed  that  the  rocket  trajectory  is 
specified  in  a  right-hand  topocentric  coordinate  system  (x,y,z). 
The  positive  x-axis  points  east,  and  the  positive  y-axis  points 
north.  The  azimuth  angle,  a,  and  the  elevation  angle,  6,  are  de¬ 
fined  in  Figure  1. 


Y  Positive  North) 


The  equations  developed  herein  are  based  upon  the  angles  9^  and  9 2 

which  are  the  components  of  the  elevation  angle  in  the  y-z  and  x-z 
planes,  respectively.  There  is,  of  course,  a  one-to-one  correspond¬ 
ence  between  the  pair  (9, a)  and  the  pair  (0p02).  The  transformation 
equations  are 


and 


9^  *  tan  ^“(tan  9  cos  a) 
92  *  tan  ^(tan  9  sin  a) 


a 

9 


tan  ^“(tan 
tan  *[tan 
tan  *(tan 


e2  cot  9^) 

9^  cot  *(tan  92  cot  0^)  ] 
9. / cos  a) . 


(1) 


(2) 


DISCUSSION 


The  launch  angles  required  to  compensate  for  wind  effect  depend  upon 
the  desired  trajectory  objective.  These  angles  are  determined  by  an 
iterative  process  of  trajectory  simulation.  The  first  estimate  is 
determined  from  the  ballistic  wind.  A  trajectory  is  computed  using 
these  angles.  If  the  simulated  objective  is  not  within  tolerance  of 
nominal,  a  correction  for  launcher  angles  is  determined  from  the  error 
vector  and  a  new  trajectory  is  simulated.  The  iteration  is  continued 
until  the  tolerance  is  achieved.  This  procedure  requires  two  algorithms, 
one  for  the  first  estimate  and  another  for  the  iteration  step. 

The  first  estimate  is  based  upon  the  ballistic  wind.  The  algorithm 
for  this  estimate  can  be  expressed  in  functional  form  by 


A61  -  fl(Wx’Wy> 
A62  -  f2(Wx>Wy) 


(3) 


where  W  and  W  are  the  East-West,  North-South  ballistic  winds,  respec¬ 
tively,  Xand  A 8^  and  A92  are  required  changes  from  nominal  settings.* 

The  functional  form  for  the  iteration  step  is  similar, 


A6^  *  (a, p) 

A©2  ■  g2(a»B) 


(4) 


where  a  and  8  are  the  component  deviations  from  the  nominal  objective 
and  A0^  and  A82  are  the  changes  in  the  launcher  settings  req  lired  for 

the  next  iteration. 


*  When  the  procedure  is  used  in  a  real-time  meteorological  system  where 
launcher  settings  are  computed  for  consecutive  (timewise)  wind  profiles, 
W^  and  Wy  are  replaced  by  the  changes  in  the  ballistic  winds  from  the 

last  profile  and  A0^  and  A02  become  changes  from  the  settings  obtained 

for  that  profile. 
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Although  the  development  presented  herein  is  directed  toward  either  of 
the  following  trajectory  objectives: 

1)  Orientation  of  velocity  vector  at  burnout,  and 

2)  The  x,y  coordinates  of  the  impact  of  either  the  afterbody  or 
a  prescribed  booster,  it  is  apparent  that  the  procedure  is  applicable 
to  many  other  constraints. 


EXPRESSIONS  FOR  THE  ALGORITHM  FUNCTIONS 


The  development  of  the  algorithms  requires  the  determination  of  suitable 
expressions  for  f^,  f£,  and  g^»  Since  there  is  little  or  no  theo¬ 
retical  indication  of  the  nature  of  these  expressions,  an  empirical  ap¬ 
proach  is  indicated.  Various  types  of  expressions  may  be  determined  by 
a  curve  fitting  procedure  and  then  tested  for  accuracy  and  suitability. 

The  empirical  approach  may,  of  course,  yield  entirely  different  function¬ 
al  forms  for  different  missile  configurations.  The  general  procedure 
is,  however,  invariant.  It  consists  of  determining  expressions  for  g^ 
and  g2,  using  these  expressions  in  the  determination  of  expressions  for 
f^  and  and  finally,  testing  the  algorithms. 

The  following  paragraphs  will  discuss  the  expressions  for  f^,  g^  and 
for  the  Athena  missile.  Other  missile  systems  will  be  discussed  in 
later  sections. 

To  determine  approximating  expressions  for  g^  and  several  trajectory 

simulations  were  performed.  Each  of  the  trajectories  was  computed  for 
no-wind  conditions;  various  changes  in  the  launch  angles  were  used  for 
the  simulations.  The  pertinent  information  (for  the  determination  of  g^ 
and  g2)  was  extracted  from  these  simulations  and  is  presented  in  Table  I. 
An  inspection  of  these  results  suggests  that  each  of  A0^*  A02^,  AX,  and 

AY  is  approximately  linear  in  each  of  the  independent  variables  and  that 
the  contributions  of  the  independent  variables  are  approximately  additive. 

This  observation  suggests  a  least-squares  fit  of  the  generic  form 


g(cl9  z2)  -  +  Be2 


(5) 


where  (e^t  t2>  represents  either  (AO-^i  A02^)  or  (AX,  AY)  as  appropriate. 
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TABLE  I 


CHANGE  IN  BURNOUT  ATTITUDE  AND  IMPACT  FOR  SPECIFIC 
CHANGES  IN  LAUNCH  ANGLES  (REGULAR  ATHENA) 


Aex  A62 

(rad) 


Burnout 


A0 


2b 


(rad) 


Impact 

AX 


(km) 


.02 

0 

.057063 

.02 

0 

.055821 

0 

-.02 

.004602 

0 

.02 

-.005573 

.01 

.01 

.028378 

.01 

0 

-.028068 

0 

.01 

-.002662 

0 

-.01 

.002420 

.01 

-.01 

-.025603 

.01 

.01 

-.030780 

.01 

-.01 

.030747 

.01 

.01 

.025773 

-.006082 

14.48 

.007035 

-15.0/ 

-.066587 

-54.77 

.065389 

51.36 

-.003154 

7.32 

.003392 

-  7.46 

.032854 

26.14 

-.033154 

-27.00 

-.030040 

-33.84 

.036510 

18.03 

-.036046 

-2C.28 

.029452 

34.04 

AY 


34.85 

■27.98 

■11.28 

13.94 

16.63 

14.82 
6.69 
5.93 

20.98 

8.01 

10.82 
23.06 
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The  departures  from  linearity  and  additivity  may  be  investigated  by 
fitting  the  form 

2  2 

g(e1#  e2)  *  ae1  +  be^^  +  ce2  +  de2  +  ec^  (6) 


Equation  (5)  will  be  referred  to  as  the  bilinear  fit  while  (6)  will  be 
called  quadratic.  Each  expression  provides  a  good  fit  to  the  data  in 
the  sense  that  the  residuals  are  small.  The  residuals  for  Eq.  (6)  are, 
as  would  be  expected,  somewhat  smaller.  Tests  revealed  that  either  of 
these  forms  provides  a  suitable  algorithm  for  the  iteration.  These 
tests  will  be  discussed  in  detail  in  a  later  section. 

A  similar  procedure  was  used  to  determine  expressions  for  f^  and 

Trajectory  simulations  were  made  with  various  values  of  ballistic  winds 

introduced  into  the  calculations.  The  iterative  procedure  was  used  to 

determine  the  launcher  angles  required  to  compensate  for  the  effect  of 

winds.  The  results  of  these  simulations  are  shown  in  Table  II.  As 

before,  the  quadratic  and  bilinear  forms  were  fit  to  the  data.  The 

independent  variables  in  this  case  were  the  components  of  ballistic 

wind,  (W  ,  W  ) .  An  examination  of  the  residuals  indicated  that  both 
x  y 

forms  provided  a  good  fit  to  the  data. 


EVALUATION  OF  THE  PROCEDURE 


Since  the  equations  developed  in  the  preceding  section  cannot  be  veri¬ 
fied  by  theoretical  techniques,  an  empirical  test  of  their  accuracy  is 
dictated.  This  section  describes  the  tests  and  presents  the  results 
therefrom. 

Sixteen  different  wind  profiles  were  used  in  the  evaluation.  Some  of 
these  profiles  were  measured  during  missile  support  operations  at  White 
Sands  Missile  Range  while  the  others  are  hypothetical.  These  wind  data 
are  shown  in  Figures  2,  3,  4  and  5.  It  should  be  observed  that  there  is 
considerable  variation  among  the  profiles;  hence,  it  is  believed  that, 
collectively,  they  provide  a  good  test  for  the  evaluation  of  the  algo¬ 
rithms. 

To  show  that  the  algorithms  are  not  restricted  to  a  given  missile  con¬ 
figuration  and  to  investigate  the  possibility  of  using  a  missile-invariant 
form  for  the  algorithm  expressions,  iour  different  rockets  were  used  in 
the  evaluation.  These  were:  two  configurations  of  the  Athena  which  have 
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TABLE 

II. 

CHANGE  IN  LAUNCH  ANGLE  REQUIRED  TO  COMPENSATE  FOR 

SPECIFIED  BALLISTIC  WIND  TO 

MATCH  EITHER  NOMINAL 

IMPACT  OR  BURNOUT  ATTITUDE 

(REGULAR  ATHENA) 

Burnout 

Impact 

w 

X 

W 

y 

A0^ 

A02 

A61  A62 

(mps) 

(rad) 

(rad) 

5 

0 

-.000210 

-.034690 

-.000702  -.034263 

-  5 

0 

.000287 

.034524 

.000404  .034279 

0 

5 

-.033055 

-.000243 

-.032452  .000455 

0 

-  5 

.033326 

. 000198 

.032233  .000621 

5 

5 

-.033193 

-.035085 

-.032940  -.034928 

5 

-  5 

.033038 

-.034346 

.031432  -.033492 

-  5 

5 

-.032857 

.034437 

-.031932  .033869 

-  5 

-  5 

.033696 

.034580 

.032871  .034731 

10 

0 

-.000350 

-.069480 

-.001356  -.068608 

-10 

0 

.000642 

.068824 

.001039  .068255 

0 

10 

-.065791 

-.000527 

-.064460  -.001021 

0 

-10 

.066858 

.000354 

.064703  .001243 
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FIGURE  2.  Wind  Profiles  Used  In  the  Evaluation. 
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FIGURE  4.  Wind  Profiles  Used  in  the  Evaluation. 
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quite  different  nominal  trajectories;  the  Aerobee-350;  and  the  Ballistic 
Missile  Target  System  (BMTS) .  A  short  description  of  each  will  be  given 
in  the  following  paragraphs. 

Two  independent  iterations  were  performed.  The  first  was  to  match  nom¬ 
inal  impact;  the  second  was  to  match  nominal  burnout  angles.  A  converg¬ 
ence  tolerance  of  1500  meters  was  selected  for  the  impact  iteration  ex¬ 
cept  for  one  case,  the  "regular"  Athena  where  a  tolerance  of  300  meters 
was  used.  This  smaller  tolerance  was  used  to  show  the  speed  with  which 
the  iteration  converges.  Compatible  tolerances,  which  differed  for  each 
rocket,  were  used  for  the  iterations  to  burnout  angles. 

The  first  and  most  exhaustive  test  was  performed  on  the  "regular"  Athena. 
This  is  a  typical  configuration  of  several  slightly  different  Athena 
missiles  fired  from  Green  River,  Utah,  to  impact  on  White  Sands  Missile 
Range.  The  Athena  is  a  two-stage  (for  the  ascent  trajectory)  unguided 
rocket  which  is  fired  at  a  nominal  elevation  angle  of  approximately 
13.5  degrees  and  achieves  an  apogee  of  approximately  250  kilometers  and 
a  range  of  725  kilometers. 

Both  the  quadratic  and  bilinear  equations  were  used  to  perform  the  itera¬ 
tions;  these  were  performed  as  independent  cases.  The  results  are  shown 
in  Tables  III  and  IV.  Two  things  should  be  observed  from  these  results. 
First,  both  expressions  provide  rapid  convergence.  Second,  there  is  no 
significant  difference  in  the  results  obtained  from  the  two  expressions. 
This  suggests  the  U3e  of  the  simpler  bilinear  equations. 

Only  the  results  from  the  use  of  the  bilinear  equations  are  presented 
for  each  of  the  other  three  test  cases.  In  every  instance  the  iteration 
converged  rapidly  -  frequently  converging  on  the  first  estimate. 

The  Aerobee-350,  used  for  the  second  test  case,  is  a  two-stage  unguided 
high-altitude  probe  which  is  fired  nearly  vertical  (elevation  angles 
between  2  and  3  degrees) .  The  nominal  range  and  altitude  vary  between 
100  and  150  km  and  320  and  400  km,  respectively.  The  results  for  this 
test  are  presented  in  Table  V. 

A  special  Athena  trajectory  requirement  was  the  subject  of  the  third 
test  case.  The  nominal  elevation  angle  for  this  trajectory  is  18.534 
degrees.  This  is  near,  but  below,  the  critical  angle  of  17.6  degrees. 
(The  critical  angle  is  defined  to  be  the  elevation  angle  which  results 
in  maximum  horizontal  range  under  no-wind  conditions.)  The  nominal 
range  for  this  trajectory  is  718.6  kilometers  while  the  maximum  (no-wind) 
range  is  721.2  kilometers.  The  nominal  impact  could  not  be  obtained  for 
wind  profiles  7,  8,  10,  and  14.  It  was  determined  that  the  maximum  range 
for  these  profiles  was  715.6,  716.4,  715.6  and  713.1  kilometers,  respec¬ 
tively.  This  decrease  in  maximum  range  is  caused  primarily  by  the 
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TABLE  III.  TEST  CASE  I,  REGULAR  ATHENA  ITERATION  TO  IMPACT 


Quadratic  Equations 


Bilinear  Equations 


Run  B.W. 

91L 

92L 

X 

Y 

61L 

82L 

X 

Y 

(mps) 

(rad) 

(km) 

(rad) 

(km) 

Nominal 

-.2121 

.1111 

306.2 

-6*43.7 

-.2121 

.1151 

306.2 

-643.7 

1  -  4.15 

-.1468 

.1479 

304.6 

-643.5 

-.1470 

.1480 

304.8 

-643.7 

-  8.93 

-.1473 

.1473 

306.2 

-643.7 

-.1472 

.1486 

306.1 

-643.7 

2  -12.91 

-.2481 

.2135 

299.9 

-647.1 

-.2473 

.2133 

299.9 

-646.0 

5.11 

-.2467 

.2155 

306.6 

-643.8 

-.2467 

.2155 

306.6 

-643.7 

-.2466 

.2153 

306.1 

-643.7 

-.2466 

.2153 

306.1 

-643.7 

3  -11.29 

-.2881 

.2011 

305.2 

-644.1 

-.2871 

.2002 

303.6 

-643.1 

10.59 

-.2879 

.2014 

306.2 

-643.7 

-.2879 

.2014 

306.2 

-643.7 

4  -12.14 

-.1602 

.2082 

300.9 

-648.2 

-.1602 

.2090 

303.0 

-647.8 

-  6.95 

-.1578 

.2095 

306.2 

-643.8 

-.1578 

.2095 

306.2 

-643.8 

5  -  1.29 

-.2274 

.1248 

307.3 

-642.5 

-.2274 

.1247 

307.3 

-642.5 

2.14 

-.2281 

.1245 

306.2 

-643.7 

-.2281 

.1245 

306.2 

-643.7 

6  2.22 

-.2075 

.0982 

306.0 

-643.4 

-.2075 

.0982 

306.0 

-643.4 

-  0.68 

7  -  2.48 

-.1731 

.1349 

299.7 

-643.3 

-.1732 

.1342 

299.8 

-643.4 

-  5.34 

-.1746 

.1378 

306.3 

-643.7 

-.1745 

.1377 

306.1 

-643.5 

8  -  3.28 

-.1569 

.1411 

304.4 

-641.7 

-.1571 

.1412 

304.5 

-641.7 

-  7.55 

-.1587 

.1423 

306.1 

-643.6 

-.1587 

.1423 

306.1 

-643.6 

9  -  5.40 

-.1703 

.1570 

306.1 

-643.7 

-.1704 

.1572 

306.6 

-643.7 

-  5.67 

-.1703 

.1570 

306.2 

-643.7 

10  -  0.13 

-.1672 

.1170 

299.5 

-644.4 

-.1672 

.1169 

299.3 

-644.5 

-  6.21 

-.1678 

.1197 

306.1 

-643.5 

-.1677 

.1196 

306.0 

-643.4 

11  0.44 

-.2272 

.1115 

306.6 

-644.7 

-.2272 

.1115 

306.6 

-644.7 

2.09 

-.2264 

.1111 

306.2 

-643.6 

-. 22S5 

.1111 

306.2 

-643.7 

12  -10.13 

-.2091 

.1925 

301.5 

-646.6 

-.2088 

.1927 

302.3 

-646.0 

-  0.26 

-.2078 

.1939 

306.3 

-643.7 

-.2078 

.1939 

306.3 

-643.7 

13  -  1.69 

-.2119 

.1281 

304.3 

-644.5 

-.2119 

.1281 

304.3 

-644.5 

0.01 

-.2116 

.1287 

306.2 

-643.7 

-.2116 

.1287 

306.2 

-643.7 

14  2.55 

-.1377 

.0972 

305.6 

-642.9 

-.1376 

.0969 

304.7 

-643.1 

-10.36 

-.1384 

.0976 

306.1 

-643.6 

-.1383 

.0976 

306.2 

-643.7 

15  -  3.66 

-.2602 

.1425 

305.8 

-644.3 

-.2600 

.1423 

305.3 

-644.2 

6.71 

-.2598 

.1426 

306.2 

-643.6 

-.2598 

.1426 

306.2 

-643.7 

16  8.10 

-.2312 

.0525 

306.5 

-641.8 

-.2311 

.0528 

307.6 

-641.5 

2.48 

-.2325 

.0528 

306.2 

-643.7 

-.2325 

.0527 

306.2 

-643.6 

13 


TABLE  IV.  TEST  CASE  I,  REGULAR  ATHENA  ITERATION  TO  BURNOUT  ANGLES  (RADIANS) 


Quadratic  Equations  Bilinear  Equations 


Run 

B.W. 

(mps) 

01L 

02L 

JD 

r-i 

CD 

02b 

01L 

02L 

6lb 

02b 

Nominal 

-.2121 

.1151 

-.6597 

.3753 

-.2121 

.1151 

-.6597 

.3753 

1 

-  4.15 

-  8.93 

-.1449 

-.1441 

.1472 

.1476 

-.6620 

-.6598 

.3745 

.3753 

-.1456 

-.1441 

.1477 

.1476 

-.6640 

-.6599 

.3762 

.3754 

2 

-12.91 

5.11 

-.2490 

-.2475 

.2145 

.2176 

-.6635 

-.6598 

.3657 

.3758 

-.2493 

-.2474 

.2147 

.2175 

-.6642 

-.6598 

.3666 

.3758 

3 

-11.29 

10.59 

-.2896 

-.2889 

.2022 

.2022 

-.6617 

-.6597 

.3755 

.3753 

-.2901 

-.2889 

.2019 

.2022 

-.6631 

-.6596 

.3748 

.3753 

4 

-12.14 
-  6.95 

-.1589 

-.1574 

.2081 

.2106 

-.6632 

-.6600 

.3674 

.3757 

-.1598 

-.1574 

.2094 

.2105 

.6663 

.6599 

.3723 

.3754 

5 

-  1.29 
2.14 

-.2279 

-.2283 

.1250 

.1244 

-.6589 

-.6598 

.3733 

.3753 

-.2279 

-.2283 

.1250 

.1244 

-.6589 

-.6598 

.3772 

.3753 

6 

2.22 
-  0.68 

-.2064 

.0976 

-.6597 

.3753 

-.2064 

.0976 

-.6597 

.3753 

7 

-  2.48 

-  5.34 

-.1721 

-.1696 

.1346 

.1361 

-.6663 

-.6601 

.3712 

.3754 

-.1723 

-.1696 

.1347 

.1361 

-.6670 

-.6602 

.3718 

.3754 

8 

-  3.28 

-  7.55 

-.1553 

-.1544 

.1406 

.1404 

-.6625 

-.6599 

.3761 

.3754 

-.1558 

-.1544 

.1409 

.1404 

-.6639 

-.6599 

.3772 

.3754 

9 

-  5.40 

-  5.67 

-.1693 

-.1690 

.1568 

.1565 

-.6608 

-.6597 

.3764 

.3753 

-.1697 

-.1690 

.1572 

.1565 

-.6619 

-.6598 

.3778 

.3753 

10 

-  0.13 

-  6.21 

-.1657 

-.1632 

.1164 

.1181 

-. 66b6 
-.6600 

.3705 

.3753 

-.1660 

-.1632 

.1164 

.1181 

-.6673 

-.6601 

.3707 

.3753 

11 

0.44 

2.09 

-.2276 

-.2283 

.1116 

.1121 

-.6576 

-.6598 

.3732 

.3754 

-.2276 

-.2283 

.1116 

.1121 

-.6577 

-.6598 

.3733 

.3754 

12 

-10.13 
-  0,26 

-.2094 

-.2089 

.1931 

.1953 

-.6604 

-.6598 

.3682 

.3756 

-.2096 

-.2089 

.1935 

.1952 

-.6611 

-.6598 

.3698 

.3755 

13 

-  1.69 
0.01 

-.2120 

-.2188 

.1282 

.1291 

-.6602 

-.6598 

.3724 

.3753 

-.2120 

-.2118 

.1282 

.1291 

-.6602 

-.6598 

.3724 

.3753 

14 

2.55 

-10.36 

->1347 

-.1332 

.0959 

.0953 

-.6641 

-.6599 

.3779 

.3754 

-.1354 

-.1333 

.0959 

.0953 

-.6658 

-.6600 

.3779 

.3754 

15 

3.66 

6.71 

-.2614 

-.2622 

.1432 

.1440 

-.6595 

-.6598 

.3723 

.3753 

-.2616 

-.2622 

.1431 

.1441 

-.6580 

-.6598 

.3720 

.3753 

16 

8.10 

2.48 

-.2307 

-.2310 

.0519 

.0513 

-.6590 

-.6598 

.3770 

.3753 

-.2310 

-.2310 

.0523 

.0513 

-.6599 

-.6597 

.3786 

.3752 

14 


TABLE  V.  APPLICATION  OF  BILINEAR  EQUATIONS  TO  AEROBEE  350 
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increased  drag  and  lift  due  to  a  head  wind.  The  results  of  this  test 
are  presented  in  Table  VI. 

The  final  test  case  was  for  the  BKTS  vehicle.  This  is  a  two-stage 
unguided  rocket  fired  at  a  nominal  elevation  angle  of  30.2  degrees 
(for  the  case  considered  here),  resulting  in  nominal  range  and  alti¬ 
tude  of  293  kilometers  and  89. A  kilometers,  respectively.  The  nominal 
elevation  angle  is  considerably  below  critical  angle.  The  results  of 
this  test  are  presented  in  Table  VII. 


CONCLUSIONS 


A  rapidly  converging  iterative  technique  to  compute  wind-compensated 
launcher  settings  for  unguided  rockets  was  developed  and  discussed. 
This  technique  is  applicable  to  a  wide  range  of  possible  trajectory 
objectives-specif ied  attitude  at  a  given  time,  space  point,  specified 
booster  impact,  etc. 

The  procedure  consists  of  two  algorithms.  The  first  provides  an  in¬ 
itial  estimate  for  the  launcher  setting  as  a  function  of  the  ballistic 
wind,  and  as  such,  it  is  easily  adaptable  to  field  operations.  The 
second  algorithm,  the  iterative  step,  converges  rapidly  (one  or  two 
passes)  to  the  desired  objective. 

Sixteen  wind  profiles  were  used  in  four  separate  test  cases  to  pro¬ 
vide  an  evaluation  of  the  procedure.  Two  different  functional  forms 
for  the  algorithms  were  investigated  -  a  quadratic  form  and  a  bilinear 
form.  Two  specific  trajectory  objectives  were  considered  in  each  test 
case  -  nominal  burnout  attitude  and  nominal  impact.  In  each  test  the 
iteration  converged  rapidly  to  the  desired  objective.  There  appeared 
to  be  no  significant  difference  between  the  results  obtained  from  the 
quadratic  and  bilinear  equations. 
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TABLE  VI.  TEST  CASE  3,  APPLICATION  OF  BILINEAR  EQUATIONS  TO  LOW  ANGLE  ATHENA 


To  Burnout  To  Impact 


Run  B.W. 

01L 

e2L 

eib 

e2b 

01L 

02L 

X 

Y 

(mps) 

(rad) 

(rad) 

(rad) 

(km) 

Nominal 

-.2871 

.1577 

-.8818 

.5732 

-.2871 

.1151 

325.1 

-640.4 

1  -  4.23 
-  8.56 

-.2190 

-.2158 

.1912 

.1923 

-.8901 

-.8823 

.5723 

.5735 

-.2063 

-.1985 

.1844 

.1822 

320.3 

324.1 

-640.7 

-639.4 

2  -11.79 
5.04 

-.3194 

-.3193 

.2574 

.2600 

-.8812 

-.8820 

.5646 

.5741 

-.3238 

-.3309 

.2586 

.2628 

325.5 

325.8 

-643.4 

-639.4 

3  -11.29 
10.77 

-.3639 

-.3607 

.2542 

.2548 

-.8901 

-.8815 

.5739 

.5733 

-.3650 

-.3693 

.2549 

.2583 

323.9 

325.2 

-643.1 

-640.2 

4  -11.14 

-  7.05 

-.2268 

-.2253 

.2495 

.2497 

-.8858 

-.8819 

.5738 

.5733 

-.2179 

-.2254 

-.2283 

.2432 

.2485 

.2492 

324.2 

326.9 

325.5 

-644.4 

-640.9 

-640.7 

5  -  1.52 
2.08 

-.3024 

.1705 

-.8852 

.5759 

-.3044 

.1716 

324.5 

-639.9 

6  2.17 

-  .43 

-.2852 

.1385 

-.8844 

.5716 

-.2835 

-.2805 

.1377 

.1370 

323.4 

324.9 

-640.7 

-640.5 

7  -  2.26 
-  4.31 

-.2527 

-.2508 

.1754 

.1818 

-.8842 

-.8820 

.5532 

.5735 

Solution  does 

not  exist 

8  -  3.54 

-  6.96 

-.2317 

-.2297 

.1856 

.1865 

-.8868 

-.8822 

.5717 

.5734 

Solution  does 

not  exist 

9  -  5.49 
-  5.69 

-.2403 

-.2384 

.2023 

.2008 

-.8877 

-.8821 

.5801 

.5731 

-.2343 

.1986 

326.0 

-640.1 

10  .16 
-  5.23 

-.2470 

-.2451 

.1549 

.1617 

-.8840 

-.8820 

.5519 

.5733 

Solution  does 

not  exist 

11  .63 

1.60 

-.2999 

-.3015 

.1523 

.1512 

-.8779 

-.8820 

.5755 

.5732 

-.3006 

-.3109 

.1527 

.1555 

331.3 

324.4 

-641.5 
-640. 4 

12  -  9.37 
-  .38 

-.2791 

-.2809 

.2359 

.2382 

-.8760 

-.8820 

.5640 

.5736 

-.2811 

-.2902 

.2358 

.2416 

325.2 

325.5 

-644.7 

-640.1 

13  -  1.42 
.01 

-.2865 

-.2867 

.1692 

.1709 

-.8806 

-.8819 

.5674 

.5733 

-.2872 

-.2908 

.1695 

.1721 

324.5 

325.2 

-642.5 

-640.2 

14  2.25 

-  9.71 

-.2138 

-.2084 

.1365 

.1385 

-.8954 

-.8827 

.5711 

.5733 

Solution  does 

not  exist 

15  -  3.36 
6.50 

-.3354 

-.3368 

.1868 

.1872 

-.8779 

-.8820 

.5707 

.5732 

-.3384 

-.3501 

.1887 

.1933 

329.7 

324.4 

-643.3 

-639.7 

16  8.01 
2.37 

-.3099 

.0904 

-.8843 

.5707 

-.3052 

.0879 

324.9 

-640.1 
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TABLE  VII.  TEST  CASE  IV,  APPLICATION  OF  BILINEAR  EQUATIONS  TO  BKTS 


To  Burnout  To  Impact 


Run 

B.W. 

eiL 

92L 

6  ib 

6  2b 

61L 

02L 

X 

Y 

(mps) 

(rad) 

(rad' 

(rad) 

(km) 

Nominal 

-.3965 

.3842 

-.6056 

.5837 

-.3965 

.3842 

199.2 

-215.1 

1 

-  6.37 

-  6.37 

-.3580 

-.3513 

.4286 

.4279 

-.6153 

-.6056 

.5859 

.5843 

-.3543 

.4292 

198.3 

-214.1 

2 

-11.12 

3.56 

-.4301 

-.4177 

.4613 

.4650 

-.6231 

-.6054 

.5808 

.5858 

-.4371 

-.4418 

.4759 

.4845 

198.6 

199.0 

-217.2 

-215.1 

3 

-12.43 

12.12 

-.4900 

-.4738 

.4701 

.4688 

-.6293 

-.6054 

.5885 

.5854 

-.5052 

-.4931 

.4939 

.4848 

197.3 

199.4 

-214.5 

-215.4 

4 

-  8.11 
-  5.41 

-.3660 

-.3572 

.4407 

.4399 

-.6184 

-.6058 

.5863 

.5846 

-.3638 

.4439 

200.1 

-216.3 

5 

-  1.88 
1.C2 

-.4051 

.3972 

-.6067 

.5853 

-.4066 

.4001 

199.7 

-214.9 

6 

0.56 
-  1.54 

-.3854 

.3804 

-.6047 

.5831 

-.3939 

.3784 

199.2 

-215.2 

7 

-  0.62 
-  1.11 

-.3894 

.3886 

-.6031 

.5826 

-.3886 

-.3810 

.3883 

.3779 

199.2 

198.7 

-213.1 

-214.5 

8 

-  3.53 

-  4.46 

-.3688 

.4089 

-.6095 

.5818 

-.3659 

.4083 

198.0 

-214.5 

9 

-  4.81 

-  5.97 

-.3594 

-.3547 

.4178 

.4154 

-.6129 

-.6058 

.5879 

.5839 

-.3556 

.4172 

200.0 

-215.4 

10 

2.17 
-  1.86 

-.3819 

-.3880 

.3692 

.3690 

-.5968 

-.6058 

.5827 

.5833 

-.3795 

-.3762 

.3652 

.3566 

200.3 

198.5 

-212.8 

-214.2 

11 

2.13 
-  1.40 

-.3851 

.3694 

-.6035 

.5837 

-.3831 

.3660 

200.5 

-216.7 

12 

-  6.63 

-  0.67 

-.3974 

-.3900 

.4303 

.4300 

-.6159 

-.6057 

.5813 

.5845 

-.3990 

-.4061 

.4364 

.4462 

199.2 

199.3 

-214.0 

-215.2 

13 

-  1.69 
0.05 

-.3965 

.3959 

-.6088 

.5821 

-.3988 

.3977 

199.1 

-216.1 

14 

-  0.25 

-  7.18 

-.3474 

.3862 

-.6053 

.5836 

-.3410 

-.3296 

.3798 

.3686 

198.1 

198.4 

-213.6 

-214.1 

15 

-  1.84 
4.18 

-.4267 

.3969 

-.6088 

.5838 

-.4311 

-.4453 

.4026 

.4161 

200.7 

199.2 

-216.9 

-215.2 

16 

6.20 

0.25 

-.3931 

-.3988 

.3411 

.3396 

-.5975 

-.6056 

.5847 

.5834 

-.3913 

.  3350 

199.4 

-214.6 

18 


BIBLIOGRAPHY 


Duncan,  Louis  D.,  and  Bernard  F.  Engebos,  “Techniques  for  Computing 
Launcher  Settings  for  Unguided  Rockets,”  ECOM-5077,  Atmospheric 
Sciences  Laboratory,  White  Sands  Missile  Range,  New  Mexico,  1966. 

Engebos,  Bernard  F. ,  and  Louis  D.  Duncan,  “A  Nomogram  for  Field  Deter¬ 
mination  of  Launcher  Angles  for  Unguided  Rockets,"  EC0M-5088, 
Atmospheric  Sciences  Laboratory,  White  Sands  Missile  Range,  New 
Mexico,  1966. 

Hennigh,  Keith  E.,  “Field  Wind  Weighting  and  Impact  Prediction  for 
Unguided  Rockets,”  NASA,  Washington,  D.  C.,  1964. 

James,  Robert  L.,  Jr.,  and  Ronald  J.  Harris,  “Calculation  of  Wind 

Compensation  for  Launching  of  Unguided  Rockets,”  Langley  Research 
Center,  Langley  Field,  Virginia,  1960. 

Lewis,  J.  V.,  “The  Effect  of  Wind  and  Rotation  of  the  Earth  on  Un¬ 
guided  Rockets,”  Ballistic  Research  Laboratory,  Aberdeen  Proving 
Ground,  Maryland,  1949. 

Walter,  Dr.  Everett  L.,  "Six-Variable  Ballistic  Model  for  a  Rocket," 
MM-445,  Missile  Meteorology  Division,  U.  S.  Army  Signal  Missile 
Support  Agency,  White  Sands  Missile  Range,  New  Mexico,  1962. 

Wilson,  G.  G.,  "A  Technique  for  Predicting  Nonlinear  Wind  Compensation 
of  Ballistic  Rocket  Systems,”  Sandia  Laboratories,  Albuquerque, 
New  Mexico,  1967. 


19 


APPENDIX  A 


BALLISTIC  WIND  WEIGHTING 


The  technique  of  wind  weighting  is  the  usual  method  of  determining 
first-order  effects  of  the  wind  on  the  trajectory  of  an  unguided 
rocket.  One  of  the  earliest  discussions  of  this  procedure  was  given 
by  Lewis  (1949)  who  also  presented  a  procedure  for  calculating  the 
information  required  by  the  weighting  technique.  Since  Lewis'  orig¬ 
inal  paper,  many  improvements  have  been  made  in  the  trajectory  simu¬ 
lation  algorithms  required  to  calculate  the  constants  for  the  wind 
weighting.  A  discussion  of  the  algorithms  for  simulation  would  re¬ 
quire  a  lengthy  discourse  of  flight  dynamics  and  aerodynamics  and 
will,  therefore,  be  omitted.  A  partial  bibliographical  listing  of 
the  theoretical  development  of  these  algorithms  is  presented  in  the 
references. 

This  appendix  will  discuss  the  techniques,  assumptions,  and  limita¬ 
tions  of  wind  weighting  and  present  the  basic  procedure  for  computing 
the  required  constants.  The  availability  of  a  computer  simulation 
program  capable  of  determining  the  trajectory  for  a  given  set  of 
conditions  is  assumed. 

Wind  velocity  is  ordinarily  a  function  of  altitude.  A  wind  profile 

is  defined  to  be  the  triple  (W  (z),  W  (z),  z)  where  z  is  altitude 

and  W  (z) ,  W  (z)  are  the  x  andXy  components  of  wind  at  height  z. 
x  y 

Wind  weighting  is  based  upon  two  fundamental  concepts  -  the  ballistic 
wind  and  the  unit  wind  effect.  The  ballistic  wind  is  that  constant 
wind  (as  a  function  of  altitude)  which  has  the  same  first-order  effect 
on  the  rocket  (trajectory)  as  the  actual  wind  profile.  The  unit  wind 
effect  is  the  wind  effect  of  a  constant  unit  wind.  A  concept  closely 
associated  with,  and  required  for  the  computation  of,  the  ballistic 
wind  is  the  wind  weighting  curve.  Let  f(z)  be  the  wind  effect  for  a 
wind  profile  consisting  of  a  uniform  unit  wind  up  to  altitude  z  and 
a  zero  wind  above  z.  The  wind  weighting  function,  g(z),  is  f(z) 
divided  by  the  unit  wind  effect;  the  graph  of  g(z)  is  the  wind  weight¬ 
ing  curve.  The  derivative  g'(z)  ■  dg(z)/dz  is  the  wind  weighting 
factor  function. 

The  ballistic  wind  is  determined  from  the  wind  profile  by  the  integral* 

BWx  *  /  W  r(z)g'  (z)dz 

BW  *  f  W  (z)g' (z)dz. 

y  *  y  ° 
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In  practice  the  integrals  are  usually  replaced  by  finite  sums.  The 

atmosphere  is  divided  into  horizontal  strata  (wind  layers)  wherein  the 

wind  is  assumed  constant;  the  average  value  over  the  layer  is  usually 

used.  If  the  boundaries  of  the  strata  are  z.  <  z0  <  ...  <  z  ,  then 

i  2  n 

the  ballistic  wind  is  approximated  by 

BWX  -  T  V«(Zi+l>  - 

i-1 

BWy  -  -  8(Zi)] 

The  value  g(z^+.p  “  gC5^)  is  often  called  the  weighting  factor  or  the 
ballistic  factor  for  the  ith  layer. 

The  wind  effect  is  calculated  by  multiplying  the  ballistic  wind  by  the 
unit  wind  effect.  These  calculations  inately  assume 

(1)  The  effect  of  the  wind  in  any  particular  stratum  is  directly 
proportional  to  the  wind  in  that  stratum. 

(2)  The  effects  of  the  wind  in  the  various  strata  are  independent 

It  is  easy  to  take  issue  with  these  assumptions;  in  fact,  one  can  pro¬ 
duce  results  of  trajectory  simulations  which  contradict  the  assumptions 
However,  for  determining  first-crder  effects,  which  is  the  aim  of  the 
wind  weighting  technique,  the  assumptions  are  reasonable. 


DETERMINATION  OF  THE  WIND  WEIGHTING  CURVE 


In  the  preceding  section  a  definition  of  wind  effect  was  omitted.  Such 
a  definition  is  unnecessary  unless  one  is  interested  in  quantitative 
results.  If  one  considers  the  trajectory  to  be  defined  by  a  family  of 
parameters,  then  wind  effect  can  be  broadly  defined  as  the  change  in 
these  parameters  due  to  the  wind  encountered  along  the  trajectory. 
Typical  parameters  are:  first-stage  impact,  second-stage  impact,  ori¬ 
entation  of  velocity  vector  at  burnout  of  a  particular  stage,  etc. 

The  calculation  of  the  wind  weighting  function  requires  the  computation 
of  a  number  of  trajectories.  Each  trajectory  simulation  is  based  upon 
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the  same  set  of  initial  conditions;  however,  different  wind  profiles  are 
used  for  each  simulation.  These  calculations  yield  a  series  of  values 
for  the  trajectory  parameters  from  which  one  determines  the  wind  weight¬ 
ing  function. 

There  are  numerous  procedures  which  can  be  followed  in  performing  the 

calculations.  The  following  procedure  is  somewhat  typical  and  serves 

as  a  more  precise  explanation  of  the  general  procedure  given  above.  For 

simplicity  of  explanation  it  will  be  assumed  that  one  is  interested  in 

wind  effect  on  impact  (either  afterbody  or  a  given  booster).  Consider 

altitudes  0  *  z  <  z,  <  z0  <  . . .  <  z  and,  for  each  0  <  i  <  n,  wind  pro- 
o  1  l  n 

files  Pi  where  each  P  consists  of  a  zero  wind  above  altitude  z i  and  a 

uniform  wind  (of  constant  nonzero  magnitude,  say  C,  and  constant  direc¬ 
tion)  below  z^.  (The  same  value  of  the  wind  velocity  is  used  in  the 
nonzero  portion  of  each  profile.)  Note  that  P  is  the  zero  profile; 
the  corresponding  trajectory  is  usually  called  the  no-wind  (nominal) 
trajectory.  Let  Ri  be  the  range,  from  launcher  to  impact,  for  the 
simulated  trajectory  which  uses  the  ith  profile.  The  unit  wind  effect 
is 


a  «  (R  -  R  )/C. 
n  o 

Specific  values  for  the  wind  weighting  function,  g(z),  at  the  points  z^ 
are 


g(Zi)  -  (R.  -  Ro)/e. 


Observe  that  g(z  )  =  0  and  g(z  )  =  1. 

o  n 

The  wind  weighting  curve  may  be  plotted  from  the  values  of  g(z^).  Of 
course,  one  must  choose  the  z.'s  sufficiently  dense  to  describe  the 
graph.  The  choice  of  z  sometimes  presents  problems.  This  altitude 
must  be  such  that  any  w?nd  effect  above  z  can  be  ignored.  Apogee 
altitude  will  clearly  satisfy  this  requirement. 


TWO-DIMENSIONAL  REFINEMENTS 


In  the  general  discussion  of  wind  weighting,  it  was  assumed  that  the 
weighting  factor  curve  and  the  unit  wind  effect  were  one-dimensional 
in  the  sense  that  the  magnitude  of  the  wind  effect  for  the  uniform 
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profiles  was  assumed  to  be  independent  of  wind  direction.  This  assump¬ 
tion  is  quite  good  for  near-vertical  firings;  however,  for  other  cases, 
the  results  can  sometimes  be  improved  by  extending  the  technique  and 
using  separate  wind  weighting  curves  and  unit  wind  effects  for  the  x 
and  y  wind  components.  The  definitions  are  similar  to  those  given 
earlier.  The  x-unit  wind  effect,  a  ,  is  the  x-component  of  the  wind 

effect  for  a  profile  with  components  W  *  1,  W  *  0  throughout.  For 

x  y 

a  given  altitude,  z,  let  P(z)  be  the  wind  profile  with  components  W  «  1, 

W  *  0  below  z  and  W  -  W  *  0  above  z.  The  wind  weighting  function  for 
y  x  y  o  s 

the  x-component,  g  (z),  is  the  wind  effect  for  P(z)  divided  by  o  .  The 

X  X 

functions  g^(z)  and  are  defined  analogously. 

The  component  wind  effects  and  wind  weighting  curves  are  calculated 
similarly  to  those  for  the  one-dimensional  case  and  will  not  be  ex¬ 
plicitly  outlined. 
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